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Abstract 

The ro-vibrational spectra of molecules dissolved in liquid 4 He nanodroplets display rotational 
structure. Where resolved, this structure has been used to determine a temperature that has been 
assumed to equal that of the intrinsic excitations of the helium droplets containing the molecules. 
Consideration of the density of states as a function of energy and total angular momentum demon- 
strates that there is a small but significant bias of the rotor populations that make the temperature 
extracted from a fit to its rotational level populations slightly higher than the temperature of the 
ripplons of the droplet. This bias grows with both the total angular momentum of the droplet and 
with the moment of inertia of the solute molecule. 
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INTRODUCTION 



A hallmark of ro-vibrational spectroscopy of molecules in liquid 4 He nanodroplets has 
been the observation of rotational structure, as in the gas phase The distribution 

of intensity of transitions in this structure has allowed the determination of the rotational 
temperature of these molecules, and it has been assumed that this provides a measure of 
the temperature of the droplets themselves This assumption has been supported 

by the close similarity of the temperatures determined for a wide range of molecules and 
by the assumption that on the long time scales between pickup of solutes and spectroscopic 
interrogation, the entire system will come to equilibrium as the droplets cool by evaporation. 

Adriaan Dockter and I have modeled the evaporative cooling of pure and doped helium 
nanodroplets using an angular momentum conserving statistical rate model [4]. These cal- 
culations have predict that the droplets cool to a temperature close to those inferred from 
experiments and previous evaporative cooling calculations, but with a vastly broader dis- 
tribution of energy and total angular momentum than for a canonical distribution at the 
same temperature. We also found, much to our initial surprise, that while the populations 
of the rotational levels of a solvated rotor follow a thermal distribution, the fitted temper- 
ature of the rotor populations is higher than that of the droplet itself. While this would 
seem to violate a fundamental principle of thermodynamics, as will be demonstrated below, 
this divergence of the two temperatures is in fact a consequence of the constraints imposed 
by conservation of angular momentum. The temperature of the helium is defined by the 
inverse of the derivative of the log of the density of states with respect to energy at fixed 
total angular momentum. However, when a rotor is excited to a state with rotor angular 
momentum quantum number j, the internal degrees of freedom of the droplet can have an- 
gular momentum values between \ J — j\ and J + j, where J is the total angular momentum 
quantum number, which is treated as a conserved quantity. Thus the derivatives of the 
density of droplet states with respect to J contributes to the rotor population. 

STATE DISTRIBUTIONS 

At the temperature of helium nanodroplets, 0.38 K, the only intrinsic droplet excitations 
are surface ripplon modes [^j]. The density of states and other thermodynamic quantities 
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can be written as a function of a reduced energy, E = E/Er where Er = 3.77 fcgK/v N 
and N is the number of helium atoms in the droplet Q|. We denote the density of states as a 
function of reduced energy and total angular momentum as pju(E, J), and the density states 
summed over angular momentum states as Pr(E) = X)j(2 J + i)p&j(E, J). Using methods 
recently described these can be well approximated by the analytical expressions 

p R (E) = exp (aE A ' 7 + bE 1 ' 7 ) (1) 

p RJ (E,J) = PR (E)(2J + l)^p^exp(-/3(E)(J+l/2) 2 ) (2) 
/3(E) = cE~ s ' 7 + dE- 13 ' 7 (3) 

if the values a = 2.5088, b = -4.3180, c = 0.8642 and d = -0.3524 are used Q. Note that 
for fixed J, the density is the number of distinct states, i.e. does not include the (2 J + 1) 
spatial quantization degeneracy of each such state. We define the ripplon microcanonical 
temperature, T&(E), and with fixed total energy and angular momentum, T R j(E, J), by 

1 =rf .u(pHW) 4 1 
T n (E) dE 7 7 K ' 



T RJ (£,J) T R (£) V 2 /W )\ dE , 

where in both cases, temperature is measured in units of reduced ripplon energy, E R , divided 
by Boltzmann's constant. 

We now consider a droplet that has a solvated rigid linear rotor, with effective rotational 
constant B, given again in units of E R . Let j be the rotational quantum number of the 
rotor, J R the rotational quantum number of the ripplons, and J the total rotational quantum 
number. The total density of states for fixed values of E, J is given by 

J+j 

p(E,J) = Y: E PRj(E-Bj(j + l),J n ) (6) 

3 Jtt=\J-j\ 

P(j), the fraction of states with rotor quantum number j, is given by 

p(E, J) 

The temperature, T(E, J), of the coupled system, at fixed E and J is given by 

1 = d\n(p(E,J)) = gj Eff|^ PRJ (g ~ gjg + 1), Jr) ^JjE) 
T{E,J) dE p{E,J) 1 ' 



A rotor temperature, Tj(E, J) can be defined from a Boltzmann fit to the rotor level pop- 
ulations P(j). This latter corresponds to the "droplet temperature" that has been reported 
in numerous experiments. This assignment of the rotor temperature to that of entire system, 
i.e. that Tj(E, J) = T(E, J) is natural assumption for a system in microcanonical equilib- 
rium, i.e. by equating P(j) to the ratio of density of states. If we ignore angular momentum 
conservation and sum over all values of J, then the above identification would be exact as 
long as the mean energy in the rotor, T r , if much less than the total energy in the droplet, E, 
so that the "heat bath" of ripplon states does not change temperature significantly over the 
range of significantly populated rotor states. Equivalently, we expect Tj(E, J) = T(E, J) if 
that Pk(E) can be approximated by Pr(E — Bj(j + 1)) = Pr(E) exp [—Bj(J + l)/T R (£ , )j. 



Note that if the finite heat capacity of the droplet heat bath is considered, then we would 
expect the population in higher rotational levels to fall off faster than predicted by the 
Boltzmann distribution, and thus the effective temperature determined by a fit to the rotor 
populations to be lower than the droplet temperature. 

When we look at the J conserving ensemble, we must consider the fact that the average 
value of Jr(Jr + 1) equals to J(J + 1) — j(j + 1) for fixed J, j. Thus the population of 
different rotor levels will be influenced not only by Tr, but also by the J dependent factor 
in prj, which has the same form as the rotational distribution function of a spherical top 
in a canonical ensemble. In order to test the size of the expected bias in the estimate of 
T(E) by Tj, we will examine numerical results for droplet of N = 3 000 and 10 000 helium 
atoms, for which Er = 0.069 and 0.0377 K respectively. Based upon our evaporative cooling 
calculations [4], we will assume that the droplets cool to a final temperature of 0.35 K, but 
have a wide range of final angular momentum values, J. This condition allows calculation 
of an isothermal curve E(J) for each droplet size. We further assume that B = 1 GHz, 
which equals 0.70 (1.27) in reduced units for droplets of 3 000 (10 000) helium atoms. This 
is the effective B constant for SFr dissolved in helium, from which the temperature of helium 

n 

nanodroplets was first measured [1|. It is found that T(E(J), J) = Trj(E(J), J) to better 
than 1 % for all values of J over the range from 0-4000. It is found that the values of 
P(j) fit a Boltzmann distribution to high accuracy, so that values of T r are well defined. 
However, over this range of J, we find a small be systematic difference in the values of 
between Tj(E(J), J) and T(E(J), J), as demonstrated in figure 1. The ratio of Tj(E(J), J) 
to T(E(J), J) is found to be approximately proportional to B~ l , i.e. is linearly proportional 
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to the effective rotational moment of inertia of the rotor in helium (including the contribution 
due to motion of the helium). As is evident from the figure, the size of the bias decreases 
for larger droplets, though less than linearly with number of helium atoms in the droplet. 

The present results demonstrate that for helium droplets with significant trapped angular 
momentum, the population of rotational levels of rotors will be biased from that predicted 
by a Boltzmann distribution at the temperature of either the droplet excitation modes or 
the temperature of the entire system. The size of the bias is on the order that should be 
detectable by careful measurement of the droplet rotational constant as the temperature of 
the pickup gas (and thus the initial collisional angular momentum) is increased. Previous 
experiments have reported a range of rotational temperatures, but up to know these do not 
appear to have been considered physically significant. For example, in a study of the cluster 
size dependence of the spectrum of SFg, the temperature was found to be constant to within 
statistical noise (0.36 — 0.40 K) for droplets of more than a few thousand He atoms, but to 
systematically increase for smaller droplets, down to a value of 0.48(20) K for droplets of a 
few hundred He atoms 0\. The authors concluded that the temperature of smaller droplets 
was higher, but the present result suggest that another possibility is that the effect they 
observed could be, at least in part, due to the bias effect predicted by the present work. 
Nauta and Miller found that for the spectrum of HF dimer, the rotational temperature was 
0.34(1) K 0, below the excepted droplet temperature of 0.38(2) K. It should be pointed out 
that the rotational constant of HF dimer in helium is three times larger than that of SFg, and 
that the pickup of two HF molecules should deposit considerably less angular momentum 
than the pickup of SFg. It should be noted that the size of the bias is a function only of the 
total rotational quantum number, and thus the effect should be essentially unchanged for a 
symmetric top or spherical top molecule. For an asymmetric top, this effect should lead to a 
failure of the transition intensities to follow that expected for a gas in thermal equilibrium if 
this bias effect is ignored. Stimulated by the present analysis, Roger Miller and collaborators 
have made an attempt to look in a systematic way at the temperature extracted from fits to 
rotor populations for different molecules. Unfortunately, the asymmetries of the lineshapes 
often observed for the lowest rovibrational transitions lead to model dependent biases in the 
fitted rotational temperatures that could not be disentangled from the expected effects [10]. 
This give yet one more reason for the need to address this problem of lineshapes of ro- 
vibrational transitions in helium nanodroples beyond the one previous attempt (U|. That 
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previous work assumed a canonical ensemble of internal states which the recent evaporative 
cooling work has now shown to be a poor assumption. 
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FIG. 1: Plot of the ratio of the rotor temperature, T r (E, J) to the droplet temperature, T(E, J) 
for a droplets of 3 000 (dashed curve) and 10 000 (solid curve) helium atoms and a solvated rotor 
with B = 1 GHz as a function of the total angular momentum quantum number, J. For each value 
of J, E has been fixed by the condition that Trj = 0.35 K 
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